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Abstract 

An analysis is given of the decay r — > fi + 7 in MSSM extensions with a vector like 
generation. Here mixing with the mirrors allows the possibility of this decay. The analysis 
is done at one loop with the exchange of charginos and neutralinos and of sleptons and 
mirror sleptons in the loops. It is shown that a branching ratio B(r — > /ry) in the range 
4.4 x 10~ 8 — 10 -9 can be gotten which would be accessible to improved experiment such as 
at SuperB factories for this decay. The effects of CP violation on this decay are also analyzed. 
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1 Introduction 



Violation of lepton flavor is an important indicator of new physics beyond the standard 
model. In the absence of a CKM type matrix in the leptonic sector, flavor violations can 
only arise due to new physics and thus decays such as U — > Ij'-f (i ^ j) are important probes 
of new physics. We focus here on the decay r — > fi + 7 on which Babar Collaboration pQ 
and Bell Collaboration [2] have put new limits on the branching ratio. Thus The current 
experimental limit on the branching ratio of this process from the BaBar Collaboration [1] 
based on 470fb~ 1 of data and from the Belle Collaboration [2] using 535 fb _1 of data is 

B{r ->■ fi + 7) < 4.4 x KT 8 at 90% CL (BaBar) 
B{r fi + 7) < 4.5 x KT 8 at 90% CL (Belle) (1) 

At the SuperB factories^ HJ E] (for a review see [6]) the limit is expected to reach B{r — > 
yU + 7) ~ 1 x 10 -9 as shown in Fig.(l). Thus it is of interest to see if theoretical estimates for 
this branching ratio lie close to the current experimental limits to be detectable in improved 
experiment. 

Here we explore this process in the presence of a new vector like generation in an extension 
of MSSM. Vector like multiplets arise quite naturally in a variety of grand unified models [7] 
and some of them can escape supermassive mass growth and can remain light down to the 
electroweak scale. Recently an analysis was given of the EDM of the tau in the framework of 
an extension of the minimal supersymmetric standard model with a vector like multiplets [8] . 
Specifically mixing of the standard model leptons with the mirror leptons, and mixing of the 
sleptons with mirror sleptons, were considered and it was found that such contributions could 
put the tau EDM in the detectable range. Here we extend this analysis to investigate the 
contributions from a vector like lepton multiplet to the flavor changing process r — > fi + 7. 
This decay is forbidden at the tree level due to vector current conservation and can only 
arise at the loop level. The current work is a logical extension of the previous works where 
mixings with a vector like multiplet and with mirrors were considered [HJ [Hll El EH 112] . 
Implications of additional vector multiplets in other contexts have been explored by many 
previous authors (see, e.g., [131 EU EH EH])- Several studies already exist on the analysis of 
r -»■ /17 decay [TTl UHl LISl EQl EH [22l |23l El]. However, none of them explore the class of 
models discussed here. 
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Figure 1: A display of the upper limits on the branching ratio B{r~ — > /j _ 7) (and for 
r~ — > /J _ 7, from the previous experiments and for the anticipated experiments as 

a function of the integrated luminosity. Figure is taken from Ref. [I]. 



2 Extension of MSSM with a Vector Multiplet 

We begin with a brief discussion on extension of MSSM where we include vector like lepton 
multiplets since such a combination is anomaly free. First under SU(3)c X SU(2)l X U(l)y 
the leptons of the three generations transform as follows 

Al= ~(l,2,-^),&~ (1,M), (1,1,0), 2 = 1,2,3 (2) 

where the last entry on the right hand side of each ~ is the value of the hypercharge Y 
defined so that Q = T3 + Y. These leptons have V — A interactions. We can now add 
a vector like multiplet where we have a fourth family of leptons with V — A interactions 
whose transformations can be gotten from Eq.(|2]) by letting i run from 1-4. A vector like 
lepton multiplet also has mirrors and so we consider these mirror leptons which have V + A 
interactions. Their quantum numbers are as follows 

(ivf) ~(1.2,^),^l~ (1,1, -l),iV L ~ (1,1,0). (3) 
The MSSM Higgs doublets as usual have the quantum numbers 

H ^{%)~^-b- H ' s {m)~ {1 - 2 - 1 2 ) - (4) 

As mentioned already we assume that the vector multiplet escapes acquiring mass at the 
GUT scale and remains light down to the electroweak scale. As in the analysis of Ref. [8] 
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interesting new physics arises when we consider the mixing of the first three generations of 
leptons with the mirrors of the vector like multiplet. Actually we will limit ourselves to the 
second and third generations since only these are relevant for the computation of the decay 
r — > /ry. Thus the superpotential of the model may be written in the form 

+h 1 Hifi llL iMl + h!^WA\ + •/:>'< A"<-1 + UaX C %L + + &Z L N L (5) 

where tpL stands for tp 3L and -ip^i stands for ty 2L . Here we assume a mixing between the 
mirror generation and the third lepton generation through the couplings f 3 , f± and f$. We 
also assume mixing between the mirror generation and the second lepton generation through 
the couplings f 3 , f' A and f' b . The above six mass terms are responsible for generating lepton 
flavor changing process. We will focus here on the supersymmetric sector. Then through the 
terms f 3 , f$, fy, f' A , f' b one can have a mixing between the third generation and the second 
generation leptons which allows the decay of r — >■ /ry through loop corrections that include 
charginos, neutralinos and scalar lepton exchanges with the photon being emitted by the 
chargino (see the left diagram of Fig. (2)) or by a charged slepton (see the right diagram of 
Fig. (2)). The mass terms for the leptons and mirrors arise from the term 

1 d 2 W 

where ip and A stand for generic two-component fermion and scalar fields. After spontaneous 
breaking of the electroweak symmetry, (< Hi >= Vi/y/2 and < H\ >= V2/V2), we have 
the following set of mass terms written in the 4-component spinor notation 



ff 1V1/V2 U 

-C m = (f R E R p, R ) f 3 f'v 2 /V2 f 3 

V f 4 Mi/A 

/ f[v 2 /y/2 h \ 

+ (v T R N R D„ R )\ -f 3 f 2Vl /V2 -f 3 \\N L \+H.c. (7) 

V f b h[v 2 /V2j ' 

Here the mass matrices are not Hermitian and one needs to use bi-unitrary transformations 
to diagonalize them. Thus we write the linear transformations 
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such that 

(hvi/V2 u o \ 

D T J \ f 3 f 2 v 2 /\/2 D T L = diag(m n ,m T!2 ,m T3 ). (9) 

V fi h lVl /V2j 

The same holds for the neutrino mass matrix 

^ /f[v 2 /V2 h \ 

D R \ ~h i^i/v 7 ^ -fl \D L = diag(m Vl ,m^,m V2 ). (10) 

V ft h[v 2 /V2j 

In Eq.Q Ti, r 2 , r 3 are the mass eigenstates and we identify the tau lepton with the eigenstate 
1, i.e., r = Ti, and we identify r 2 with a heavy mirror eigenstate with a mass in the hundreds 
of GeV and r 3 is identified as the muon. Similarly V\, z/ 2 , z/ 3 are the mass eigenstates for the 
neutrinos, where we identify v x as the light tau neutrino, v 2 as the heavier mass eigen state 
and ^3 as the muon neutrino. 

Next we consider the mixings of the charged sleptons and the charged mirror sleptons. 
The mass 2 matrix of the slepton - mirror slepton comes from three sources, the F term, the 
D term of the potential and soft susy breaking terms. Using the superpotential of Eq.([5]) the 
mass terms arising from it after the breaking of the electroweak symmetry are given by C E 
and Co 

-C F = (ml + \f 3 \ 2 + \f,\ 2 )E R E R + (m 2 N + \f 3 \ 2 + l^ 2 )^* 
+« + \U\ 2 + \f' 4 \ 2 )E L El + (m 2 N + \f 5 \ 2 + \n\ 2 )N L Nl 
+ (m 2 + |/ 4 | 2 )f^ + (ml + \M 2 )u tR D* tR + (m\ + \h\ 2 )r L r L 

+ ( m l + I/IR/Wr + (™l + UsH/Wtz, + ( m l T + IM^tlKl 
+ K, + I/sH^l + (rnl M + |/s| 2 )^^ 
+{— m T fi* tan (3f L f R — m^fj,* tan (3N L N R — m Ur jj* cot (iv TL v* rR 

-m M // tan f3jx L jj* R - m u ^* cot /Su^l^r 
-m E fi* cot PE L E* R + (m E ft + m T f 4 )E L f* L 
+ (m E fa + m T fl)E R f* R + (m E fz + m^E^ 
+ (m E f 4 + m^f^ERp,^ + (m v J 5 - m N f^)N L i)* L 
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+(rn N f 5 - m v j;)N R v* R + {m v J' 5 - m N f*)N L i>* L 

+fsf;^Ll>rt + fctf^RKR + H.C.}. (11) 

Similarly the mass terms arising from the D term are given by 

- C D = Kn 2 z cos 2 6 W cos2(3{u rL i)* L - f L f* L + z>^z>* L - fi L p,* L + E R E* R - N R N* R } 

+ ^m| sin 2 9 W cos 2f3{u rL i>* L + f L f* L + z> m l^l + p>Lp>* L 

-E R E* R - N R N* R + 2E L E* L - 2r R r* R - 2jl R jl R }. (12) 

In addition we have the following set of soft breaking terms 

v soft = m 2 tL ^;J 1 tL + M 2 x r*r + m*J;*J; l + m^Kl 

+K D ^; L + M 2 T f?f c L + M 2 /^l + M 2 E E* L E L + M 2 N N* L N L 

+f 2 A N H[x CJ N L - f 2 A E Hlx C3 E L + H.c.} (13) 

From Cf,d and by giving the neutral Higgs their vacuum expectation values in V so f t we can 
produce the the mass 2 matrix M~ in the basis (tl, El,t r , E r , fl R ). We label the matrix 
elements of these as (M 2 )^ = M 2 - where 

M 2 U = M 2 L + m 2 r + \f 3 \ 2 - m 2 z cos2f5{ 1 - - sin 2 6 W ), 

M| 2 = M 2 + m% + |/ 4 | 2 + \f' 4 \ 2 + m 2 z cos2/3 sin 2 6 W , 
M 3 2 3 = M 2 + m 2 + |/ 4 1 2 - m 2 z cos2/3 sin 2 9 W , 
M 2 4 = M 2 + m 2 E + \f 3 \ 2 + \f' 3 \ 2 + m 2 z cos20(± - sin 2 6 W ), 

M 2 5 = M 2 L + ml+ \f 3 \ 2 - m 2 z cos2P{ 1 - - sin 2 6 W ), 
M 2 m = Ml + ml + |/ 4 | 2 - m 2 z cos2(3 sin 2 9 W , 
M 2 2 = M 2 * = m E f;+m T f A , 
M 2 3 = Mil = m T (A* T - //tan/3), 
M 2 4 = M 2 * = 0, M 2 5 = M 2 * = / 3 /*, 
M 2 * = M 6 2 * = 0, M 2 3 = M 3 2 * = 0, 
M 2 4 = M42 = m E (A* E - /i cot 0), M 2 5 = M52 = m E f 3 + mj',*, 



5 



M| 6 = Ml; = 0, Ml = Ml* = m E f 4 + m T f*, M 2 5 = M 2 * = 0, M 3 2 6 = M 6 2 * = f 4 f; 

Ml = Mil = 0, M 2 6 = Mg4 = m E f: + m^, 

M 2 6 = Ml; = m,(A; - Man/3) (14) 

Here the terms M 2 1; M 2 3 , M| 1; M| 3 arise from soft breaking in the sector ft, f^, the terms 
M| 5 , M| 6 , Mg 5 , Mg 6 arise from soft breaking in the sector and the terms M 2 2 ,M 24 , 

M 2 2 ,M% 4 arise from soft breaking in the sector El,Er. The other terms arise from mix- 
ing between the staus, smuons and the mirrors. We assume that all the masses are of the 
electroweak size so all the terms enter in the mass 2 matrix. We diagonalize this hermitian 
mass 2 matrix by the unitary transformation D T "<M?D T = diag(M? v M? , M? , M?, M?, M? ). 
There is a similar mass 2 matrix in the sneutrino sector. In the basis [y T L-, A"l, z> tR , N r , z/^, z> mR ) 
we can write the sneutrino mass 2 matrix in the form (M~)ij = m?- where 

m 2 u = M 2 rL + m 2 UT + \f 3 \ 2 + ^m 2 z cos2f3, 
m\ 2 = M 2 N + m 2 N + |/ 5 | 2 + |/ 5 | 2 , ml, = M 2 Ut + m 2 T + |/ 5 | 2 , 
m 2 4 = M 2 + m 2 N + \f 3 \ 2 + |^| 2 - l -m 2 z cos2f3, 

ml, = M 2 L + m 2 M + Iftf + l -m 2 z cos2P, 

<6 = Ml+ml + \ftf, 

m i2 = m l* = m u T h - "W 3 *> 
m i3 = m ii = m vXK T - fJ> cot 0), m 2 14 = m 2 4 \ = 0, 

2 2* n 2 2* ri r* 2 2* n 

m 14 = m 41 = 0, m 15 = m 51 = / 3 / 3 , m 16 = m 61 = 0, 

m 23 = m 32 = ° ? m 24 = m 42 = ™> N (A* N - fl tan ft) , 77l| 5 = TO52 = -771jv / 3 + m^/g* , 

m 26 = m 62 = °> m 34 = m 43 = m N fb ~ ™vJt, 

2 2* n 2 2* r ri* 2 2* n 

m 35 = m 53 = °> m 36 = m 63 = /s/s » m 45 = m 54 = 

m 2 6 = m 2 * = -m^/g + mjv/s*, m 2 6 = m 2 * = m V)i {A*^ - //cot /3). (15) 

As in the charged slepton sector here also the terms m 2 1; m 2 3 , m 31 , m 33 arise from soft break- 
ing in the sector i> r £, ^vr, the terms mf 5 , m 2 6 , mg 5 , arise from soft breaking in the sector 
i / nL,v flR , and the terms m 22 ,m 24 , m 42> m 44 arise from soft breaking in the sector Nl,Nr. 
The other terms arise from mixing between the physical sector and the mirror sector. Again 
as in the charged lepton sector we assume that all the masses are of the electroweak size 
so all the terms enter in the mass 2 matrix. This mass 2 matrix can be diagonalized by the 
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unitary transformation D U ^M?D U = diag(M? i} M? 2 , M? 3 , Mf 4 , M? , M|J. The physical tau 
and neutrino states are r = Ti, v t = u±, and the states r 2 , ^2 are heavy states with mostly 
mirror particle content. Similarly /1 = T3, = z/3. The states fj, 2 = 1 — 6 are the slepton 
and sneutrino mass eigenstates. 



where 



3 Interactions of Charginos and Neutralinos 

The chargino exchange contribution to the decay of the tau into a muon and a photon arises 
through the left loop diagram of Fig. (2). The relevant part of Lagrangian that generates 
this contribution is given by 

3 2 6 

- £ r _,_ x + = E E E f «[ c &d p L + CaijPR&j + H.c. (16) 

a=l i=l j=l 

+V ll Dl\ a D\ J + - KeVaDZfi^l (17) 

where D u is the diagonalizing matrix of the scalar 6x6 mass 2 matrix for the scalar neutrino 
as defined above. k n , n T etc that enter in the equation above are defined by 

{K N , Kr, «„) = — =— -, (KB, Ku) = -/=— — ■ (18) 

V 2Mw cos [3 \/2M w sm p 



In Eq.(17) U and V are the matrices that diagonalize the chargino mass matrix Mc so that 

WMcV- 1 = dtag(m±,m±). (19) 



The neutralino exchange contribution to the tau decay arises through the right loop 
diagram of Fig. (2). The relevant part of Lagrangian that generates this contribution is 
given by 

3 4 6 

- £ T _ f _ x o = E E E UC%Pt + C'* 3 P r \ X \t 3 + H.c, (20) 

a=l i=l j=l 

7 



T ' 
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Figure 2: The diagrams that allow decay of the r into \i + 7 via supersymmetric loops 
involving the chargino and the sneutrino (left) and the neutralino and the stau (right) with 
emission of the photon from the charged particle inside the loop. 



where as stated earlier r = t% and /i = T3. In Eq.(20) C' L and C' R are defined by 



where D T is the diagonalizing matrix of the 6x6 slepton mass 2 matrix. 

gm E X* 4] , g , 1 

"*> ~ 2m w sinp" 0E > " + ( 2 

,„ #sin 2 w , gm E X Aj 

— eA y — a 2 -, o s — — ; — -, 

J cos t>w 2mw sin p 



IE, 



X 2j (^- -sin 2 w ), 



and 

gm T X 3j 

a • = — 

TJ 2mw cos (3 



1 



gm^X 3j . 



(21) 



(22) 



" x;*(-^ + sin 2 ^), 

,(23) 



where 



osin 2 ^ grrirX^ gm^X*. 

7ri 7w- eA lj+ cog ^ a 2j ., o Tj - 2mwCOg/3 > <W- 2m w ,cos/3 

= (Xij cos 9 W + X 2 j sin 8 W ), X' 2j = (—Xij sin 9 W + X 2 j cos 0vp), 



(24) 
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and where the matrix X diagonlizes the neutralino mass matrix so that 

X T M x oX = diag(m x o i ,m x o 2 ,m x o 3 ,m x o 4 ). 



(25) 



4 The analysis of r — > ji + 7 Decay Width 

The decay r — > is induced by one-loop electric and magnetic transition dipole moments, 
which arise from the diagrams of Fig. (2). In the dipole moment loop, the incoming muon 
is replaced by a tau lepton. For an incoming tau of momentum p and a resulting muon of 
momentum p', we define the amplitude 

< fi(j/)\J a \T(p) >= U^p')V a U T {p) (26) 

where 

rjg) = F l^^l + m^ml + (27) 

with q = p' — p and where rrif denotes the mass of the fermion /. The branching ratio of 
r — > fi + 7 is given by 

= 5C | m ;K + %)- " f7(0|p + K ' m (28) 

where the form factors and arise from the chargino and the neutralino contributions 
as follows 

F^(0) = F% + + 

F^(0) = F^ + +F^ o (29) 
The chargino contribution F^ + is given by 

2 6 / v yiy-2 

jrTix _ \^ y^ { rn T {m T + m^) {r<L R R 

*2x+ ~ 1^ U 647r 2 m 2 ^^lij + G 3ii<-lijKll^^J 

1=1 j = l Xi + Xi + 

+ (30) 



where 



= —^——{-2x 3 -3x 2 + 6x-l + 6x 2 lnx} (31) 

OyCC L J 
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and 

F 2 (x) = - — l —r^{Zx 2 - 4x + 1 - 2x 2 lnx} (32) 
(x — 1) 

The neutralino contribution F^ is given by 



4 6 _ ( _i_ ^ M 2 
V "2^2J 1927r2m 2 ' ( ! ^ 3 fe 

i=l j=l Xi Xi 



. m T (m T + m M ) /L /L+ n iR n iR* 



(m T + rrifj) ( /L >R# r ,iR r ,,L*^ F / \i / Q q\ 
" 64vr%>~o K :1 " - + ^ }F4( ^ )J (33) 

Xi 

where 

F 3 (x) = -, --{-x 3 + Qx 2 - 3x - 2 - Qxlnx} (34) 

(x — 1) 

and 

F 4 (x)= ^ {-x 2 + 1 + 2xlnx} (35) 
(x 1) 

The chargino contribution is given by 

2 6/ , \ TfP' 

pT/i \^ \^ + 'm^)m iii + L m R , x -.+ 

^3 X + - 2^ 2^ 32vr 2 m? i G 3ii G i« ~ C 3 l j C iij^5{j^) (3b) 

i=l j=l v i v i 

where 

F »<*>=2(^l)i<-* + 3+ T^> < 37 » 
The neutralino contribution F^ Q is given by 

4 6/ . \ ™ 2 

- 2^ 2^ 327r 2 M 2_ l C 3ij G l« - C Kj C Uj 1^6(^2") (38) 

i=l j=l T i T i 

where 

*^^T? {x + l + — } (39) 



5 Estimate of size of B(r — > fij) 

In this section we give a numerical analysis of B(t — > /ry) for the model where we include a 
leptonic vector multiplet. As discussed in the previous sections the flavor changing processes 
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Figure 3: An exhibition of the depen- 
dence of B(t — > ^7) on m when uin = 
120, m E = 150, |/ 3 | = \&\ =90, |/ 4 | = 
\f A \ =100, |M = \&\ =80, \A Q \ =150, 
mi = 50, m 2 = 100, /x = 150, % 3 = 
Xs =0-6, X4 = X4 =0-4, Xs = x' 5 =0-6, 
as =0.5, ajy =0.8, and tan/3 =5, 10, 15, 
20 (in ascending order at mo = 300). Here 
and in Figs. (4-7) masses are in GeV and 
angles are in rad. 



Figure 4: An exhibition of the de- 
pendence of B(t — > /ry) on |/ 3 | when 
m =900, rriN = 150, m E = 180, 

\ f >\ =ioo, \u\ = \fi\ =ioo, \f 5 \ = 

\f[.\ =70, \A \ =100, mi = 50, m 2 = 100, 
A* = 150, x 3 = x's =0-6, X4 = X4 =°- 4 > 
Xs = X5 =0.6, a E =0.5, «at =0.8, and 
tan/3 =5, 10, 15, 20 (in ascending order 
at I/3 1 = 100.) 



arise from the mixings between the standard model leptons and the mirrors in the vector mul- 
tiplet. The mixing matrices between leptons and mirrors are diagonalized using bi-unitary 
transformations with matrices D T R and D T L . The input parameters for this sector of the pa- 
rameter space are m T , m E , m^, / 3 , / 4 , / 3 , where / 3 , f±, / 3 and f A are complex masses with 
CP violating phases x 3 , X4> X3 Xa- F° r the slepton mass 2 matrices we need the extra input 
parameters of the susy breaking sector, Mr^, M E , M T , M x , M^ L , M M , A T ,A E , A^,A N ,^, tan /3. 
For the sneutrino mass 2 matrices we have more input parameters, Mn, M„ r , M u , A Vfi , An, A Ue 
m N-, /sj /s- The chargino and neutralino sectors need the extra two parameters mi,m 2 . In 
the analysis we will include phases since dipole moments are sensitive to phases (for a review 
see [25J). Here for simplicity we assume that the only parameters that are complex in the 
above matrix elements are A E , An, A t , A^, A v , / 5 and f 5 which have the phases a E , a at, 
a T , a M , a v , Xs and Xs- To simplify the analysis we set a u = = a T = 0. Thus the CP 
violating phases that would play a role in this analysis are 

X3,X4,X5,X3iXi,X5i a E,aZN. (40) 

With the above in mind, the electric dipole moments of the electron, the neutron and of 
the Hg atom vanish and we do not need to worry about them satisfying their upper limit 
constraints. To reduce the number of input parameters we assume equality of the scalar 
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Figure 5: An exhibition of the de- 
pendence of B(t — > /ry) on x 3 when 
tan/3 =10, m N = 170, m E = 200, 
|/ 3 | = l/al =250, |/ 4 | = \fi\ =400, |/ 5 | = 
\f' 5 \ =90, \A \ =130, rhi = 90, m 2 = 
80, /i = 120, x'3 =0-8, X4 = Xa =0-9, 
Xs = X5 =1-6, «e =1.0, Oat =0.9, and 
m =900, 800, 700, 600, 500 (in ascend- 
ing order at % 3 = 0.0.) 



Figure 6: An exhibition of the de- 
pendence of B(t — > fij) on X4 when 
m =800, tan/3 =15, m^ = 160, ttie = 
220, \&\ =150, |/ 4 | = \fi\ =200, |/ 5 | = 
\f^\ =100, |A | =160, mi = 100, m 2 = 
90, fi = 150, xs = Xs =0-6, X4 =0-8, 
X5 = X5 =1-0, oe =-4, aw =0.8, and 
l/s I =300, 250, 200, 150 (in ascending or- 
der at Xi = 0.0.) 



masses and of the trilinear couplings so that M a = m , a = t l , E, t, x, v, a*, /j-l, N and \ A { \ = 
\A \, i = E,N,r,u,fi. 

Fig. (3) gives an analysis of B(r — > /ry) as a function of mo for values of tan /3 = 5, 10, 15, 20 
with other inputs as given in the caption of Fig. (3). The branching ratio depends on the 
chargino and neutralino exchange contributions to F 2 and F 3 defined in Eq.(29 ) which depend 
on m through the slepton masses that enter the loops. Fig. (3) exhibits a sharp dependence 
on tan /3 which enters F2 and F3 also via the slepton masses as well as through the chargino 
and neutralino mass matrices. Further, the couplings C L,R and C' L,R are also affected by 
variations in m and tan/3. The analysis of Fig. (3) shows that there is a significant part of 
the parameter space where Bir — > /ry) lies in the range O(10~ 8 ) consistent with the upper 
limit of Eq.g. Fig. (4) gives an analysis of Bir — > ^7) as a function of I/3I, where f 3 is an 
off diagonal term in the mass matrix of Eq.Q, for tan (3 values as in Fig. ^ and the other 
inputs are as given in the caption of Fig. (4). As in Fig. (3) one finds a sharp dependence 
on tan/3. This dependence of I/3I arises since it enters in the matrix elements diagonalizing 
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Figure 7: An exhibition of the dependence of B(r — > /ry) on tan (3 when mo =700, = 200, 
m E = 300, |/ 3 | = \f,\ =180, |/ 4 | = =100, |/ 5 | = \fl\ =150, \A Q \ =360, m x = 120, 
m 2 = 80, fi = 140, X3 = X 3 =0 - 7 > X4 = X4 =0.9, Xs = X5 =- 6 > «s =-9, ajv =0.4, and 
X3 =1.2, 0.8, 0.5, 0.1 (in ascending order at tan/3 = 30.) 



matrices D T L R and this way it affects both chargino and neutralino exchange contributions. 
The entire parameter space exhibited in this figure is consistent with the upper limits of 
Eq.Q. 

We discuss now the effect of CP phases on B{r — > fij). As mentioned above the phases 
of Eq.(40) have no effect on the EDMs of the electron, on the EDM of the neutron or on 
EDM of the Hg atom and these phases only affect phenomena related to the second and 
the third generation leptons. Fig. (5) gives a display of B(t — > /ry) as a function of X3 f° r 
values of mo = 900, 800, 700, 600, 500 GeV (in ascending order) when tan (3 = 10 and the 
other inputs are as shown in the caption of Fig. (5). Here one finds that Bir — > /ry) has a 
significant dependence on % 3 . Thus, for instance, for the case m = 500 GeV (top curve) 
one finds that B(t — > /ry) can vary in the range (1 x 10~ 8 — 4 x 10~ 8 ) as X3 varies in the 
range (0, 7r). Again B{r — > /ry) displayed in this analysis is consistent with the upper limit 
of Eq.([T]) over the entire range of parameters exhibited. 

Another analysis on the dependence of Bir — > fi'j) on CP phases is exhibited in Fig. (6) 
where a plot of B{t — > yiry) as a function of %4 is given for the case when |/ 3 | = (300,250, 
200, 150) GeV (in ascending order), tan = 15 and other inputs are as given in the caption 
of Fig. (6). Again a very significant variation in Bir — > /ry) is seen as \i varies in the 
range (0,7r). Specifically one finds that for the case I/3I = 150, B{r — > yLry) varies in the 
range (8 x 10 -9 — 3 x 10~ 8 ). Further, over the entire parameter space analysed in Fig. (6) 
B{t — > fij) is consistent with the upper limit of Eq.Q. Finally, in Fig. (7) we exhibit the 
dependence of B{r — > ^7) on tan/3 when xz = 1-2, 0.8, 0.5, 0.1 (in ascending order) with other 
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parameters as defined in the caption of Fig. (7). A sharp dependence of B(t —> /ry) on tan/3 
can be seen. Specifically one finds that for the case X3 = 0.1 (the top curve) B(r —> fij) varies 
in the range (1 x 1CT 10 — 3 x 1CT 8 ) which is more than an order of magnitude variation as 
tan/3 varies in the range (5-30). 

In summary in the analyses presented in Fig. (3-7), one finds that B(t — > yrf) can be quite 
large often lying just below the current experimental limits which implies that this part of 
the parameter space will be accessible to future experiments, specifically SuperB factories 
which can probe B(r — > /I'y) as low as 10~ 9 . We note that the flavor changing interactions 
of Eq.([5]) also contribute to the muon anomalous magnetic moment g M — 2 which is very 
precisely determined experimentally. This can come about by the exchange of a tau and a 
photon in the loop but since each vertex is one loop order, the contribution is three loop 
order which would be tiny compared to other standard model electroweak contributions. 

6 Conclusion 

Lepton flavor changing processes provide an important window to new physics beyond the 
standard model. In this work we have analyzed the decay r — > fi + 7 in extensions of the 
MSSM with vector like leptonic multiplets which are anomaly free. Specifically we con- 
sider mixings between the standard model generations of leptons with the mirror leptons in 
the vector multiplet. It is because of these mixings which are parametrized by ^3, /4, /s and 
f'z-i f'i-, f§ as defined in Eq. |5p that lepton flavor violations appear. We focus on the supersym- 
metric sector and compute contributions to this process arising from diagrams with exchange 
of charginos and sneutrinos in the loop and with the exchange of neutralinos and staus in the 
loop. These loops do not preserve lepton flavor. A full analytic analysis of these loops was 
given which constitute the main result of this work. A numerical analysis was also carried 
out and it is found that there exists a significant part of the parameter space where one can 
have the branching ratio for this process in the range 4.4 x 10~ 8 — 10~ 9 , where 4.4 x 1CT 8 
at 90% CL is the upper limit from BaBar (see Eq. ([!])) and the lower limit is the sensitivity 
that the SuperB factories will achieve. Thus it is very likely that improved experiment with 
a better sensitivity may be able to probe this class of models. 
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